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Abstract 

Equilibrium relationships involving solids are based in bulk thermodynamic properties that concern ideal crystals of infinite size. 
However, real processes towards equilibrium imply nucleation and development of finite molecular- scale entities. The configuration 
of these early-stage clusters and the estimation of their excess energies with respect to the ideal crystal are key to understand the 
macroscopic behavior of a given system. Here, starting from the ideal atomic positions in bulk calcite, aragonite, or vaterite, the 
relaxation of finite clusters of CaC03 in vacuum has been explored. With the aim of determining the impact of the cluster size on 
its energy and its geometrical configuration, a series of CaC03 clusters have been simulated and their lattice energies calculated. 
The cluster geometry has been fully optimized at constant pressure and its energy has been determined using GULP. A wide variety 
of clusters ranging from 1 to 2000 formulae has been considered for each starting structure (calcite, aragonite, or vaterite). GULP 
calculations have been carried out using the pair potentials set derived by Rohl etal^. In a number of cases the final configuration 
has been checked using the DFT code SIESTA ElO, and an excellent agreement has been found. Although these simulations do 
not represent fully realistic scenarios (U|5l|6]|, some results are relevant from the point of view of the polymorphic precipitation 
of CaC03 . A thorough analysis of the diffraction patterns of the relaxed clusters has allowed addressing the fundamental question 
about the critical size that a cluster should have to be considered as truly calcite, aragonite, or vaterite. 
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1. Introduction 

Equilibrium relationships involving solids are based on bulk 
thermodynamic properties that concern ideal crystals of infinite 
size. However, real processes towards equilibrium imply devel- 
opment of finite molecular-scale entities. The configuration of 
these early-stage clusters and the estimation of their excess en- 
ergies with respect to the ideal crystal are keys to understand- 
ing the macroscopic behaviour of a given system. As nucle- 
ation events are difficult to study experimentally, both because 
they occur spontaneously and because the nucleus size is very 
small, atomistic simulations are a suitable tool for understand- 
ing the early stages of crystallisation. Here, starting from the 
ideal atomic positions in calcite and aragonite, the relaxation in 
vacuum of finite clusters of CaC03 is explored. Nucleation and 
growth of calcium carbonate phases constitute a very impor- 
tant subject of research in a wide variety of fields. A complete 
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study of CaCC>3 should include many different aspects: size 
and shape of the critical nuclei under diverse conditions, possi- 
bility of nucleation from precursor phases, nucleus energy and 
nucleus surface energy, relationship nucleus- substrate in het- 
erogeneous nucleation We present a preliminary study of nu- 
cleation of calcium carbonate where nuclei are considered to be 
isolated from any previous phase or substrate. Even when this 
situation does in no way represent realistic conditions, it can 
be a helpful first approach to more complex studies. In ambi- 
ent conditions of pressure and temperature, calcite is the most 
stable phase of CaC03. However, if the Ostwald-Lussac law 
of phases were strictly true, the crystallisation of calcite would 
begin with nucleation of an amorphous material, which would 
then transform into vaterite, then aragonite, and finally calcite. 
Actually, numerous experimental works have shown the exis- 
tence of precursor phases on calcite crystallisation. 

2. Calculations: method and procedure 

With the aim of determining the influence of the cluster size 
on its energy and on its geometrical configuration, a series of 
virtual CaC03 clusters was built and their lattice energy calcu- 
lated. For each simulation run the starting configuration was 
an isolated spherical and electrically neutral cluster of calcium 
carbonate (calcite or aragonite) containing a definite number of 
CaC03 formulae. The cluster geometry was fully optimized at 
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Table 1: Parameters of the force field. Carbonate anions are defined as molecules, within which Coulomb interactions are excluded and C/O atoms are deemed to be 
bonded. All atoms are represented by a core, except O atoms, which also have a shell. The charges are: Ca core = +2.0, C core = +1.343539, O core = +1.018487, 
and O shell = -2.133. 



constant pressure and its energy was determined using the GULP 
program (dill). A wide variety of clusters ranging from 1 to 
2000 formulae was considered for both calcite and aragonite 
initial configurations For these calculations atomistic simula- 
tion techniques based on the Born model of solids were em- 
ployed. In this method, the atoms are considered to be charged 
balls, which are free to interact with each other. These interac- 
tions are long-range electrostatic forces and short-range forces 
which can be described using simple analytical functions. The 
components of the short-range forces include both the repul- 
sions and the van der Waals attractions between neighbouring 
electron charge clouds. GULP calculations were carried out us- 
ing the pair potentials set derived by Rohl et al Q]]. This set 
of potentials includes the polarizability of the carbonate ions in 
the system via the shell model of Dick and Overhauser (9). The 
oxygen ions are described by a shell of zero mass representing 
the electronic charge cloud connected to a core containing all 
the mass of the ion. The total charge of the ion is the sum of 
the charges of core and shell. The position of the core repre- 
sents the position of the ion in the crystal lattice. The shell and 
core are connected by a harmonic spring with a force constant. 
All intramolecular potentials are defined to act between cores 
while intermolecular interactions act on shells, where applica- 
ble. Furthermore, all Coulomb interactions are excluded within 
molecules. Although the shell model is relatively simple, it has 
been very successful in modelling a variety of properties. This 
set of potentials is transferable between the different CaC03 
polymorphs and it was shown to be excellent at simulating sur- 
faces and environments with low coordination (l8l fT0l[TTl[T2ln 
.The parameters of the force field are given in [I] 

In order to ensure the validity of this force field in the sim- 
ulation of very small clusters, the starting configurations with 
10 or less formulae unit were also optimised and their ener- 
gies calculated with the DFT-based code SIESTA (El El). For 
these calculations the pseudopotentials for Ca, C and O by 
Archer (2004) were used and a double Z-polarized basis set was 
considered 1 13 ]. The energies and final relaxed configurations 
were verified to be consistent with the GULP results. The final 
atomic coordinates of the relaxed clusters were compared with 
the crystal structure of calcite or aragonite. With this aim, the 
X-ray diffraction patterns of all the final relaxed clusters were 



calculated with the Debye formula and compared with the cor- 
responding aragonite or calcite pattern. 
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Figure 2: Example of four different configuration of aragonite clusters with 14 
formulae. All of them have the same radius and chemical content, but their 
energies are significantly different. 



3. Results 

The calculated total energies (in kJ/molformula) of calcite 
and aragonite un-relaxed clusters from 1 to 2000 CaCC>3 for- 
mulae are represented in |l(a)| In some cases (mainly for small 
clusters with less than 30 formula units), there are diverse ini- 
tial possible configurations with the same number of formulae 
and similar radius, whose calculated energies are considerably 
different. For example, three different aragonite clusters with 
14 formula units are represented in [2] As can be seen, their 
calculated energies per mol are significantly different. In such 
cases, the energies of every possible configuration were calcu- 



lated, and the one with lower energy was plotted in 1(a) The 
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Figure 1: Energy of un-relaxed |(a)| and relaxed |(b)| clusters. The horizontal dashed lines in |(b)| represent the lattice energies of the corresponding bulks, calcite 
and aragonite. The logarithmic scale chosen for the X axis allows to observe the big differences of energy in small clusters. 



general trend is similar in both calcite and aragonite cases: there 
is a decrease in energy per formula with the increasing number 
of formulae in clusters. This increase is most important in clus- 
ters with a small number of CaC03 formulae. As is shown in 
the plots, data can roughly be fitted to logarithmic functions of 
the number of formulae, but the energy decrease is not con- 
tinuous and there are specific cluster sizes that are particularly 
stable. In these cases the increment of formulae (in a small in- 
terval) means an increase in the cluster energy. For example, 
the un-relaxed calcite cluster with 24 formulae is more stable 
than clusters with 30 and 34 formulae; or the un-relaxed arag- 
onite cluster with 10 formulae is more stable than clusters with 
12, 14 and 16 formulae. These deviations from the general ten- 
dency are most important in clusters with a small number of 
formulae. 

The calculated total energy of the relaxed clusters (in 



kJ/molformula) is represented in 1(b) The horizontal dashed 
lines below the curves represent the energies of calcite and 
aragonite bulks calculated with the same set of potentials. As 
in un-relaxed clusters, a decrease is energy per formula with 
the increasing number of formulae is observed. Here, the di- 
vergence of the data from the general tendency is significantly 
lower than in the case of un-relaxed clusters. Explored clusters 
do not reach energies as low as calcite or aragonite bulks. A 
2000 formulae relaxed cluster with calcite starting configura- 
tion is about 25 kJ/molmore energetic than the corresponding 
calcite crystal, and a 2000 formulae relaxed cluster with arag- 
onite starting configuration is about 30 kJ/molmore energetic 
than the aragonite bulk. 

4. Discussion 

For clarity in discussion, clusters will be classified in three 
groups. The first group will include small clusters, with less 
than 12 formulae. In this case, all the structural units (Ca or 
CO3) are in the cluster surface and all the coordination poly- 
hedra of Ca are incomplete. Clusters with intermediate size 
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Figure 3: Energy of small relaxed clusters (a). Example of relaxed small cluster 
with 8 formulae (b). 



between 14 and 160 CaCC>3 formulae will be considered in a 
second group. The un-relaxed clusters belonging to this second 
group have a considerable number of structural units in their 
surface, but some coordination polyhedra of Ca are complete. 
The third group will include clusters from 170 to 2000 CaCC>3 
formulae. In un-relaxed clusters of this group most of Ca com- 
plete their coordination polyhedron. 

4.1. Small clusters 

The energies of small un-relaxed and relaxed CaCC>3 clus- 
ters whose initial structure is calcite or aragonite are compared 
in[5J\. Small clusters of calcite and aragonite relax to the same 
or to equivalent relaxed configurations with the same energy. 
As a representative example, a relaxed cluster with 8 CaCC>3 
is reresented in [5J3. No order between calcium and carbonate 
can be found in relaxed small clusters. As expected, the sim- 
ulated X-ray diffraction patterns of small clusters do not show 
any significant reflection. These patterns are included in [6] 

4.2. Intermediate size clusters 

The energies of intermediate relaxed CaCC>3 clusters are 
compared in |4]\. Energies and configurations of relaxed inter- 
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Figure 4: Energy of intermediate size relaxed clusters (a). Example of relaxed 
intermediate cluster with 50 formulae (b) 




Figure 5: Energy of large relaxed clusters (a). Example of relaxed intermediate 
cluster with 300 formulae (b) 

mediate clusters are clearly different depending on their initial 
configuration. In some cases clusters relaxed from the arago- 
nite configuration are more stable and in other cases clusters 
derived from the calcite configuration are less energetic. A rep- 
resentative example of relaxed intermediate cluster is shown in 
[4j3. Relaxed intermediate clusters of calcite and aragonite are 
partially ordered. The rows of alternating Ca and CO3, charac- 
teristic of calcite and aragonite structures, can be observed in 
the inner part of the relaxed clusters. Moreover, in this inner 
part, the relative orientation of the CO3 groups is almost en- 
tirely preserved and in general, the Ca coordination polyhedra 
can still be identified although they are distorted. However, the 
structural units in the cluster surface relax to clearly different 
positions. Carbonate triangular groups tend to be positioned 
tangent to the cluster surface, and their O atoms try to com- 
plete the coordination polyhedra of the superficial Ca atoms. 
Despite the relative disorder in relaxed clusters, the simulated 
X-ray diffraction patters show peaks that can be identified as 
important reflections of the initial structure calcite or aragonite. 
Simulated X-ray diffraction patterns of intermediate clusters are 
included in [6] 

4.3. Large clusters 

The energies of large relaxed clusters of CaCC>3 are shown 
in[5j\. In large clusters, the final energy of relaxed configura- 
tions depends on the initial configuration. As can be seen, all 
relaxed clusters explored in this range have lower energy if the 
initial configuration is calcite-like. An example of a large re- 
laxed cluster is shown in [5j3. Relaxed clusters belonging to this 
group show the calcite or aragonite structure of the correspond- 
ing initial un-relaxed cluster. Except in their surface, the inter- 



atomic distances and coordination polyhedra are essentially the 
same as those of calcite or aragonite bulks. As in intermedi- 
ate clusters, the main effect on relaxation affects to the planar 
CO3 groups in the cluster surfaces, which move and rotate to be 
positioned tangent to the surface. As can be seen in[6j the simu- 
lated X-ray diffraction patterns of large clusters are very similar 
to the corresponding aragonite or calcite diffractograms. 

5. Conclusions 

Although the simulations carried out in this work do not rep- 
resent any fully realistic scenarios, we can find some relevant 
conclusions from the point of view of the polymorphic precip- 
itation of CaCC>3. Small stable clusters of calcium carbonate 
that can be interpreted as nuclei in the first stages of nucleation, 
do not present calcite or aragonite structure. For CaCC>3 entities 
below 12 formulae, the amorphous configuration is more stable 
than any other ordered phase considered. Even though calcite is 
the more stable polymorph of calcium carbonate at room tem- 
perature, some aragonite clusters in the range of 12-160 CaCC>3 
formulae are more stable tan the corresponding calcite clusters. 
One of the reasons of the existence of precursor amorphous or 
aragonite phases in the crystallization of calcite can be linked 
with these conclusions. Nevertheless, this study shows that the 
classical theory of nucleation is insufficient when an atomistic 
approach to nucleation is required. 
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Figure 6: Simulated diffraction patterns of representative caclite-like (a) and 
aragonite-like relaxed clusters (b). The main reflections of aragonite or calcite 
are shown as bars 
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